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We present asymptotically exact solutions of an incommensurate Harper equation — one-dimensional 
Schrodinger equation of one particle on a lattice in a cosine potential. The wave functions can be 
written as an infinite product of string polynomials. The roots of these polynomials are solutions 
of Bethe equations. They are classified according to the string hypothesis. The string hypothesis 
gives asymptotically exact values of roots and reveals the hierarchical structure of the spectrum of 
the Harper equation. 
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1. Incommensurate quantum systems present a broad 
variety of complicated fractal objects and strange sets. 
The Harper equation (or almost Mathieu equation) is 
one of the most notable: 



ipn+i + ipn-i + 2Acos(0 + 2^ni])vb n = Eij) r , 



(1) 



It has two competing periods r\~ x and 1, and describes 
a particle on a one-dimensional lattice in a periodic po- 
tential. The same equation (|l|) also describes a particle 
on a two-dimensional square lattice in a uniform mag- 
netic field 77 flux quanta per elementary plaquette (Azbel- 
Hofstadter problem). Then A is an anisotropy of the 
square lattice and parameter 9 is a Bloch multiplier. 

When r\ is an irrational number, i.e., the period of 
the potential is incommensurate with the period of the 
lattice, the spectrum and wave functions are known to 
have a rich hierarchical structure @- @- The physics of 
the problem depends on the value of parameter A and 
falls into three universality classes. For A < 1 all wave 
functions are extended and the energy spectrum has a 
non-zero Lebesgue measure 4(1 — A). For A > 1 all wave 
functions are localized (in this letter we consider only 
"typical" irrational numbers which can not be approxi- 
mated by rationals too well) and the energy spectrum has 
zero measure. At the most interesting "critical" point 
A = 1 wave functions neither localized nor extended, but 
exhibit a power law scaling, while the spectrum is sin- 
gular continuous (a set without isolated points but with 
zero Lebesgue measure). 

Since the empirical observations of Hofstadter j|] , the 
evidence has been mounting that these spectra are reg- 
ular and universal rather than erratic or "chaotic" . Few 
years ago it has been shown (7^[|] that despite the com- 
plexity of the spectrum, the Harper equation at any ratio- 
nal r\ — P/Q is integrable and can be "solved" by means 
of the Bethe Ansatz (BA). This had opened the possi- 
bility of describing the complex behavior of an incom- 
mensurate system as a limit of a sequence of integrable 



models. In this letter we analyze the Bethe equations |uj 
for A = 1 and show how the (BA) reveals the hierarchi- 
cal structure of the spectrum. We show that (i) solutions 
of the BA equations (roots) fall into complexes — strings; 
(ii) the number of roots in a string (length of the string), 
known as Takahashi-Suzuki numbers [TT| are intimately 
related to the Hall conductance of a state; (iii) strings 
provide an asymptotically exact wave functions for the 
entire spectrum. 

To introduce the scaling concept let us consider an ir- 
rational 77 as the limit of rational approximants 77, = -5- . 
With rational number r\j instead of 77 the potential in (fy) 
becomes commensurate with the lattice and the problem 
is reduced to diagonalization of a Qj x Qj matrix. The 
spectrum consists of Qj bands separated by Qj — 1 gaps 
p"2[ . In this letter we show that the wave function of the 
problem with rational rjj can be written approximately 
as a product of a wave function for the ?7j-i (ancestor 
wave function) and some polynomial which roots form 
a so-called string. By recursion the wave function with 
irrational rj appears as an infinite product of string poly- 
nomials. The multiplicative structure thus obtained is 
asymptotically exact (see eq. 11). It resembles the multi- 
plicative formulae conjectured in [Jl3|. 

2. The Bethe Ansatz. If rj = P/Q is a rational number 
then the solution of (|l|) can be chosen as Bloch function 
ip n = e tk * n Xn where \n = Xn+Q- Although the BA solu- 
tion is known for arbitrary values of Bloch multipliers k x 
and k y = 9 Q , below we consider only the points of Bril- 
louin zone (k x , k y ) which correspond to edges and centers 
of bands. At these (rational) points the BA solution is 
particularly simple iQ. BA solution can be obtained in 
two steps. The first step is the isospectral transforma- 
tion: 



M 



j>n = c nm ^(~ipq n 



(2) 



m=0 

to the difference equation: 
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pnzE^(z) — iq(z + irnq 2 )(z — iKq 2 )ty(qz) 

- iq^ 1 (z - iTnqi)(z + inq^)^!(zq^ 1 ), (3) 

where p = exp(z fcx+fc ^ 7rP ), M = 2Q - 1 if P is odd 
and M = Q — 1 if P is even; g = e I7r?? . The choice of 
t = ±1 yields either levels at the center or at the edges 
of bands respectively. The rational points are degenerate 
and are labeled by discrete parameters n, p = ±1. The 
coefficients in (|J) are given by "quantum dilogarithms" 

Cam- H\e q l_,W+l/2 )■ ,4! 
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The second step is based on the integrability of eq.||): 
Q its (regular) solutions are polynomials of z of the de- 
gree Q - 1: 



(5) 
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The roots {z{\ satisfy the system of Bethe equations: 
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Solutions of the BA equations give the wave functions of 
the Harper equation at band's edges and centers. Their 
energy is given by 



ipq Q (q-q x ) 
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(7) 



At first glance, the BA equations (^) look even more 
complicated than the original Harper equation. However, 
the BA equations (|^) provide a better description of the 
problem in the most interesting, incommensurate, limit 
Q — > oo. In a sense, the incommensurate limit plays the 
role of the thermodynamic limit, with Q being the size 
of the system. 

3. String hypothesis. Below we formulate the string 
hypothesis which allows us to obtain the approximate 
solutions of the BA equations. The hypothesis is based on 
the analysis of singularities of the BA, and is supported 
by extensive numerics. We present a detailed analysis of 
singularities in a more extended publication [l5| . Here 
we just formulate the string hypothesis and present some 
immediate consequences. First, we need the notion of 
strings, a Hall conductance, a hierarchical spectral flow 
and a hierarchical tree. 

a) A string of spin I, parity v = ±1 and center a; is a 
set of 21 + 1 complex numbers: 



vx{q l 2 ,q l - ,...,q l '} 



(8) 



where x is a real positive number and cfi 



21 + 1 



= ±1. 



b) The Hall conductance a(k) of the band k — 1, . . . ,Q 
is the Chern class of the band, i.e., the number of zeros of 
the wave function as a function of k x , k y in the Brillouin 
zone. We will also need the Hall conductance of k filled 
bands: Ufe. The latter is assigned to the fc-th gap, so 
that a(k) = o~k — Cfe-i- They are determined by the 
Diophantine equation Pq,p7[| 



P<Jk = k (mod Q) 



(9) 



restricted to the range — Q/2 < <Jk < Q/2. 

c) Hierarchical spectral flow. We refer to the band k' 
of the problem with parameter rf = P' /Q', as to the par- 
ent band of the parent generation corresponding to the 
daughter band k of the daughter generation rj = P/Q if: 
(i) the Hall conductance of the band k is the ratio be- 
tween the difference of the number of states of parent and 
daughter bands and the difference between the periods of 
generations 
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(ii) rj' and ry are closest rationals with Q' < Q, and 
(k' - 1)/Q' < (k - l/2)/Q < k'/Q'. These conditions 
may be viewed as the integrated Streda formula (see Ref. 
Jl8|), which determines the spectral flow while varying 
the period from 77' to rj. Starting from some particular 
band k of the generation 77 one can obtain from Eq. ( |To| ) 
the sequence of "ancestor" bands belonging to genera- 
tions « = 5 — subsequent rational approximants to 77. 
The bands belonging to different periods r\j together with 
parent-daughter relations the hierarchical tree. Let us no- 
tice that the sequence r\j differs from the approximants 
obtained by truncating the continued fraction expansion 
of rj. Instead, it consists of intermediate fractions p9| . 
Let us also notice, that parents of states of a given gen- 
eration do not necessarily belong to the same generation. 

Now we are ready to formulate the string hypothesis. 
At large Q each solution of the BA consists of strings, 
so that each state can be labeled by spins {lj, lj-i, ■ ■ .} 
and parities {vj, Vj—i, . . .} of strings, such that the total 
number of roots J2j=i(2^+1) = Q~ 1- We refer to the set 
of lengths and parities of strings constituting the solution 
for a given energy level as to a string content of this 
level. The length of the longest string in a string content 
of a given energy level is the Hall conductance of the 
corresponding energy band: 2Z+1 = |er(fc)|. The "period" 
of this string qi = e™ 71 ' is uniquely determined by the 
requirement that 77; = is the best approximant for 
the period 77 of the state. The parity of the longest string 
is vi = —iq\ +1 ^ 2 K = (— l)fa']/c. The remaining roots are 
given by the solution of the BA equation for the parent 
state of the parent generation. In other words 
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^daugther^j „ ^parcnt^ J[ (z - X^qD- (11) 
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The accuracy of eg. (|1 1|) is determined by the accuracy 
of approximation of r\ by rj and r\i which is of the or- 
der of 1/Z 2 (rj, rj' and eta; are all rational approximants 
of some "typical" irrational number). Thus the formula 
( pi] ) is asymptotically exact. It gives the roots of the 
largest string of the length 21 + 1 with the accuracy 
l/l 2 which gets higher along the hierarchical tree when 
ij — > irrational. 

The complete string composition of a given state can 
be obtained by means of the iterative procedure. Starting 
from a state of the generation 77, first we find the branch 
of the hierarchical tree that this state belong to. Then 
we determine the Hall conductances of all states of the 
branch down to the origin. This determines the lengths, 
periods and parities of strings. 

The strings hierarchy has been obtained through the 
analysis of singularities of BA (see [^5)). As it was 
expected a set of possible lengths of strings is a set 
of Takahashi- Suzuki numbers (dimensions of irreducible 
representations of U q (SL 2 ) with definite parity [pof). 
Eq.([)J) provides a relation between them and Hall con- 
ductances. 

The iterative procedure provides an algorithm for writ- 
ing down the wave functions of all states for any rational 
rj. The only unknowns are the centers of strings. They, 
however approach 1 with accuracy 0(1/ 1). The accuracy 
of the recursive eq.(ll) is 0(l~ 2 ). The string content of a 
state (i.e., lengths and parities of strings) is a topological 
characteristics, while the centers of strings are not. 

To illustrate the iterative procedure let us consider the 
bottom edges of the lowest band of the spectrum and 



choose rj 
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to be the golden mean. The sequence 



of rational approximants is given by ratios of subsequent 
Fibonacci numbers iji = — , where the Fi are Fibonacci 
numbers (F { = F^ 2 + F^i and F Q = F\ = 1). The 
set of Hall conductances = Takahashi- Suzuki numbers = 
allowed lengths of strings are again Fibonacci numbers: 
Qk = Fk- The considered branch of hierarchical tree 
connects edges (r = — l,/t=l,/z=l) of the lowest bands 
of generations 773ft = F§k-i/ Fsk- Their string content 
consists of pairs of strings with lengths 2l n + 1 = F3 n +i, 
n = 0, 1, . . . , k — 1, parities vi n = +1 and inverted centers 
xi n and x~[ l ■ According to the string hypothesis the wave 
function of this state is 



fc-i i„ 



*rX)- 



(12) 



7i=0j=-J„ 



Centers of the strings xi n are close to 1 but can not be 
obtained from the string hypothesis only. We find em- 
pirically (see Fig|l|) that for this branch they can be fit 
by xi n 1 + with £ = 1.14. Equation @ in the 

limit k — > 00 gives the wave function of the lowest energy 
level of the Harper equation with rj — (y/5— 1)/2 — golden 
mean. More accurately, it gives an asymptotically exact 
value for " /( ^' ? f +3JV at k -> 00. 




2 x 1 -string + 

2 x 5-string x 

2x21 -string * 

34/55 ansatz □ 



FIG. 1. Shown in the figure are the roots of the lowest 
rational state of the Harper equation at rj — 34/55. The 
roots are organized in strings: two of the length 1(+), two 
of the length 5(x) and two of the length 21(*). These roots 
can be approximated by ansatz (^) and are plotted in the 
Figure as boxes(|"~|). Roots tend to repel each other. This 
"splitting" is of the order of 0(1/1). 

As an another example, let us consider a branch of 
the hierarchical tree for the golden mean, which goes 
through levels r 2 k = 2F2k ~^ 2k - 1 -f 5+ ^ 1Q 1 - ) of generations 
f]2k = F 2 k-i/F 2 k- Here fe = 1, 2, . . . and r 2k is the num- 
ber of the level counted from the bottom of the spec- 
trum. This branch leads to the point of the spectrum 
characterized by an irrational number (filling factor) 
C = linifc^oo r 2 k /F 2 k = (5 — v5)/10. Each state along 
the branch consists of strings of lengths 2s n + 1 — F 2n+ i 
and parities v Sn = (— 1) 3 W 3 1, where n = 0, 1, . . . , A; — 1 
and {•} denotes fractional part. The wave function along 
the branch is 



fc-i 



*(*\v2k) « n n (* 



(13) 



n=0 j- 



The iterative procedure allows one to generate the 
string content and wave functions out of two irrational 
numbers: a period rj and a filling factor Q. 

4. Numerical evidence for the string hypothesis. We 
now present a fragment of an extensive numerical anal- 
ysis that we have performed in order to check the string 
hypothesis. Let us consider the bottom edge of the low- 
est energy band for 779 = 34/55 from the golden mean se- 
quence. According to (|l^) the string content of this level 
consists of strings of lengths (1, 1, 5, 5, 21, 21). All string 
parities in this example are +1. The actual numerically 
calculated roots are shown in Fig.J^. One notices imme- 
diately that the roots organize themselves into distinct 
groups whose members have roughly the same radius: 
two groups with 1 root(+), two with 5 roots(x), and 
two groups with 21 roots (*). The centers of strings are 
x ln = 2.10, 1.23, 0.81, 1.05, w 1 + 1.14/F^+i. The roots 
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produced by eg. (|T^) (|~| on Fig.|l|) are in a good agree- 
ment with the numerical results. In fact one can prove 
that their accuracy of each root of the string of the length 
21 + 1 is 0{l/l 2 ). 

The overlap between the asymptotic ansatz wave func- 
tion (12) and the direct numerical solution of the eq.(||) 
for ?7i5 = ||P is (^direct | ^ansatz) « 0.99475. To see how 
the accuracy of ansatz (|l2|) is increasing with the length 
of the string we replace all roots of ansatz wave func- 
tion ^ansatz except for the ones belonging to two largest 
strings of the length 377 by numerically obtained roots. 
The new wave function ^j n term is in some sense interme- 
diate between pure ansatz V'ansatz and numerical ^direct 
wave functions. The overlap (^direct IV'interm) ~ 0.999435. 
This increase of an accuracy of approximation is due to 
a higher accuracy of the long strings of the length 377 
with respect to the shorter ones. 

5. Gaps. A direct application of the string hypothesis 
is the calculation of the gap distribution p(D), i.e., the 
number of gaps with magnitude between D and D + dD. 
To find the gap distribution we estimate the size of the 
smallest gap in the spectrum of Harper equation with 
rj = P/Q. Let us, first, consider an arbitrary gap of a 
generation with P-even and trace the genealogy of lower 
and upper edges (r = —1) of the gap along two paths 
J_ and J + of the tree, until we find the nearest common 
ancestor at the generation with denominator Q' . The 
minimal gap would correspond to the shortest paths J± , 
i.e., when Q is the denominator of the generation which 
is parent of P/Q. The string decompositions of the two 
edges have a common part consisting of strings of length 
smaller than Q' (they correspond to the common path 
from the origin of the tree to the parent generation Q') 
and different strings with lengths bigger than Q' belong- 
ing to the paths J±. Then according to ^7§, the width 
of the gap is the difference between energies of strings 
with length greater than Q' along the path J_ and the 
energies of strings along J + . The energy of a string with 
a spin I which connects the center of the band (t = +1) 
with the edge of the band (r = — 1) is 



1 , 
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For a typical n the length of the largest string 21 + 1 
scales as Q along the sequence of rational approximants, 
which gives an estimate e; ~ 1/Q 2 - The difference of 
energies of strings (i.e., the width of the minimal gap) is 
of the same order. This gives the size of the minimal gap 

Anin ~ l/Q 2 . 

While moving along the tree new gaps appear, while 
older gaps stay approximately the same. There is no 
overlap between bands — all gaps are open (except for 
rj having even denominators and E = where bands 
touch) Therefore the full set of gaps in the spec- 

trum for irrational 7yj rr consists of Q — 1 gaps which has 



been opened at the period rj = P/Q plus the smaller gaps 
which open at rational approximants to ry; rr with denomi- 
nators higher than Q. One can write J D ,q-, p(D)dD = 
Q—l ~ Q. Assuming smooth power law gap distribution 
p{D) ~ D~ 7 one obtains 7 = —3/2 and the gap distribu- 
tion function: p{D) ~ D~ 3 / 2 . This confirms numerical 
analysis of Ref . [Ell . 

6. Stating the problem: scaling hypothesis. The string 
hypothesis solves the Bethe Ansatz equations with an ac- 
curacy 0{Q~ 2 ) for roots belonging to the largest strings. 
This gives asymptotically exact results for the "evolu- 
tion" along the hierarchical tree in the incommensurate 
limit Q — > 00 (see cq.(|ll|)). However, the most interest- 
ing quantitative characteristics of the spectrum are actu- 
ally in the finite size corrections of the order of Q~ 2 to 
the bare value of strings. Among them are the anomalous 
dimensions of the spectrum. Let us choose a branch J of 
the hierarchical tree and consider energies of the states 
along the branch Ej{J). The scaling hypothesis states 
that they converge to a point of the spectrum E{J) such 
that Q 2 ~^\E.j{J)-E{J)\ is bounded but does not con- 
verge to zero. The numbers e j are anomalous exponents. 
They depend on the branch (signature of multifractality) 
and on rj (according to ref. for rj — v/ ^~ 1 exponents 
ej vary between 0.171 and —0.374). Can anomalous di- 
mensions be found analytically? They are determined 
by the finite size corrections to the strings which, we be- 
lieve, can be obtained via a more detailed study of the 
BA equations. This is a technically involved but a fasci- 
nating problem. Its solution may suggest the conformal 
field theory approach, which has been proven to be ef- 
fective for finding the finite size corrections of integrable 
systems, without the actual solving the Bethe Ansatz. 
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